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Abstract The question how the spin of the nucleon is distributed among its quark and gluon constituents 
is still a subject of intense investigations. Lattice QCD has progressed to provide information about spin 
fractions and orbital angular momentum contributions for up- and down-quarks in the proton, at a typical 
scale ij,^ ~ 4GeV^. On the other hand, chiral quark models have traditionally been used for orientation 
at low momentum scales. In the comparison of such model calculations with experiment or lattice QCD, 
fixing the model scale and the treatment of scale evolution are essential. In this paper, we present a 
refined model calculation and a QCD evolution of lattice results up to next-to-next-to-leading order. We 
compare this approach with the My hrer- Thomas scenario for resolving the proton spin puzzle. 
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1 Introduction 

How is the total spin 1/2 of the nucleon distributed among 
its quark and gluon constituents? This question has been 
intensely discussed ever since the EMC experiment pre- 
sented first results for the spin asymmetry in polarized 
muon proton scattering in 1987 IL . This measurement in- 
dicated that only about 15% or less of the nucleon spin is 
built up by quark spins, although with sizeable statistical 
and systematic uncertainties. More recent measurements 
of HERMES and COMPASS |2,3 and their QCD analy- 
sis |4jT6i showed that the nucleon receives still only about 
one third of its spin from quark spins: 



Zi^H ERMES 

(5GeV^) = 
0.330 ±0.011th,„. ± 0.025, 



± 0.028, 



(1) 



determined at a scale /i^ — 5GeV^. This is in stark con- 
trast to naive model calculations, as for example in the 
non-relativistic quark model that suggests AS — 1. Rel- 
ativistic effects reduce AS to about two thirds, still far 
too large in comparison with Eq. ([Tl . Myhrer and Thomas 
proposed in (TUiOj that AS coula be further reduced by 
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including pion cloud contributions and corrections from 
one gluon exchanges. With such corrections they end up 
with a result for AS that is consistent with experiment. 
The missing ^ 60 — 70% of the nucleon spin reappear 
entirely as orbital angular momentum of up and down 
quarks. On the other hand, L^+d, appears to be in strong 
contrast to lattice calculations (lT-13 where the orbital 
angular momentum contribution Lu+d comes out close to 
13 . To explain this difference, Thomas fs] proposed 



to consider the renormalization scale (/i-) dependence of 
the quantities appearing in the nucleon spin sum rule |14| 



]^AS + Lq + Lg + AG=]^, 



(2) 



defined by the following expectation values taken for a 
spin-up state of the proton, |-P-f): 



AS 
AG 

L„ 



{P+\J £x^^^^^i^\P+), 
{P+\J d^x{E^A^ - E^A')\P+), 
{P+\i J (f'x^lj\x^d'^ -x^d^)ijj\P+), 
{P+\j d^xE'ix^d^ ~ x^d^)A'\P+). 



(3) 
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Here is the quark field, i?* and are the gluon elec- 
tric field and gauge potential. A sum over quark flavors 
is implicit in the definition of the flavor singlet quanti- 
ties in Eq. (Isl . Contributions in the non-singlet sector will 
be denoted Dy AEu-d, L^-d etc. . The Lg is the orbital 
angular momentum contribution from gluons and AG is 
the gluon spin part. It is important to note that Lg, Lg 
and AG in Eq. (|3| are not explicitly gauge invariant. A 
manifestly gauge invariant decomposition and its relation 
to moments of generalized parton distributions was pre- 



sented by Ji in 16 



-AS + L'^' + Jg 



GI 



1 

2' 



(4) 



where AS is given as before, L'^^ is obtained from Lq 
replacing 9^ by the gauge-covariant derivative, 9^ — >■ D^, 
and the total gluon angular momentum is defined as 



{P+\J d^x[xx {ExB)]3\P^ 



(5) 



Using a leading order QCD evolution of the spin contri- 
butions from the low, hadronic model scale to the higher 
scale of the lattice results, it was shown in Ref. |8 that 
it is possible to find at least a qualitative agreement with 
the lattice data. 

With these previous achievements in mind, the pur- 
pose of the present work is twofold: first, we extend the 
QCD evolution to next-to-leading (NLO) and next-to-next- 
to-leading (NNLO) order and perform a backwards evo- 
lution starting from lattice results. This approach has the 
advantage that the scale dependence of the spin contribu- 
tions is rather weak at the higher scale of lattice results, 
and that the extrapolation therefore does not suffer from 
the uncertainty of the slope in /z at low scales. Most im- 
portantly, proceeding in this way we do not have to fix the 
model scale a priori, which is generically difficult, but have 
the possibility to compare model results over a wider range 
of low scales with the downward-evolved lattice data. As 
a further extension, we use not only the perturbative cou- 
pling cks(/i) in the evolution equations but employ also 
a frequently suggested "non-perturbative" strong coupling 
that approaches a constant a^^ax ^^e infrared region. 

The second purpose of this work is to reexamine the 
model calculations of ff^TU' and also to study possible 
improvements (Section [3|. Given these results and the 
evolved lattice data, we conclude with a discussion in Sec- 
tion S] 



2 QCD evolution of lattice results 

In this section our aim is to evolve results from lattice 
QCD, usually provided in the MS scheme at a scale fi^ ~ 
4 GeV^ , down to the low scales characteristic of model cal- 
culations. The lattice calculations were performed on the 
basis of manifestly gauge invariant operators. The com- 
putations correspond to the spin decomposition proposed 
by Ji, Eq. Q. For the remainder of this section, we will 



therefore employ the gauge invariant definitions of the spin 
observables. We drop the superscript GI in the following 
for better readability. To obtain the complete set of evo- 
lution equations for all individual parts of the spin sum 
rule, we define the orbital angular momentum of quarks 
as Lq = Jq — ^AS, of gluons as Lg — Jg — AG (for dis- 
cussions of the latter definition. 



16 



see Refs. 

Note that the gauge invariant AG cannot 
sented in terms of a local operator 



18]). 

e repre- 
191, but it can be de- 



fined as the lowest x-moment of the gauge invariant gluon 
spin distribution, Ag{x). Despite remarkable experimen- 
tal and theoretical efforts with respect to polarized parton 
distributions (4||6| [20|]22] , little is known so far about the 
magnitude of AG. Concerning the numerical evaluation 
of the evolution equations, we will therefore concentrate 
on the quark spin, the quark orbital angular momentum 
and the total angular momentum of the gluons. As will be 
shown below, this can be done without explicit knowledge 
about AG and Lg = Jg — AG. It then follows that the evo- 
lution of all quantities of interest can also be performed at 
NNLO, employing known results for the relevant anoma- 
lous dimensions from the literature. 

The total angular momentum contributions Jg and Jg 
are introduced as in [16 in the framework of generalized 
parton distributions. We observe that Jg and Jg mix in 
exactly the same way under renormalization as the (sym- 
metric and traceless) quark and gluon energy momentum 
tensors. This can be seen for example by rewriting 



(p,s|j;,|p,s) = 

-e'i'' lim 
dA^^ 



dA3 



(P- 



2 " 



^,S\T^^g\P-^,S) 



i2nf6iA), (6) 



where P is the quark momentum, s the quark helicity, and 
Z\ is a momentum difference between incoming and outgo- 
ing quark. Here, the additional derivative with respect to 
the momentum transfer, A^^ , cannot have any influence on 
the singular behavior of the operators. Therefore they mix 
in the same manner. The QCD evolution equations for Jq 
and Jg are constructed using the spin-2 singlet anomalous 
dimension given at next-to-leading order in [23 24 and at 
next-to-next-to-leading order in 25 ■ 27 . This yields 



Jo 



(471) 



-0-2 



as 1 


32 
9 


" 47r 1^ 


32 
9 


— biup 




+ hiiiF 


diUp 


- b27^F - 




+ b2np - 






'd2nF + e2n^p \ ( Jq 
doHp - eoni \ J„ 



(7) 

for np flavoursfcompare also f28l), with entries ai,bi,... 
given m Table 1 For the non-singiet combination J^^ , 
we find Eq. ([8|. The evolution equations for the spin con- 
tributions at NNLO in the MS (MS) scheme (29l[30] (these 
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Table 1. Coefficients entering the evolution equations ([7|. 

ai bi di 02 &2 C2 d2 62 

11744 416 611 5514208 , 1280 ^/q\ 134888 , 2560 1136 670871 _ 2600 f-(q\ 8830 

243 81 81 6561 ~'~ 81 ' 2187 27 ^^'^t 243 4374 27 ^y^) 729 



dlnu2 " 47r 9 " \AtvJ V 243 81 



a, 1^5514208 1280 C(3) _ 167200 tif _ 1280nFC(3) _ 896 n% \ 



- 



47r/ V 6561 81 2187 27 729 



two schemes are simply connected through a change in the 



renormalization scale) are given by 31-33 

/ Snp \ / AU' 



.AnJ y^m + lriF -pi) \AG 



.47r 



3 / 472 



16ni 





735, 



(9) 



At NNLO, the anomalous dimensions and ^gg are still 
unknown, while the upper row (7^^, 7gg) has been obtained 
as described in fSSl. Here the QCD beta functions are 

/3o = ll-^, Pi^lQ2~^nF. (10) 

We emphasize that in the chosen renormalization scheme, 
the evolution of AS is independent of AG, even at NNLO. 
Furthermore, since Jq + Jg = 1/2 at any scale, the evo- 
lution of Jq does not require an independent knowledge 
of the value of Jg. Hence one finds the remarkable re- 
sult, already mentioned above, that neither AG nor Lg — 
Jg — AG are actually required in practice for the scale 
evolution of Lq = Jq~ AS /2. As a consequence the evo- 
lution of all the quantities in Eq. Q can be performed at 
NNLO. 

Employing the definitions of Lq and Lg given above, 
fully consistent coupled evolution equations for the orbital 
angular momenta of quarks and of gluons can be written. 




Table 2. Lattice QCD results from Ref. jssj for AS 12 and L, 
in the MS scheme at /^^ = 4 GeV'^, separated into u- and d- 
quark contributions. Statistical and estimated systematic un- 
certainties due to the renormalization are given in the form 

(. . .)stat(. . .) ren- 





AE/2 


L, 


u 


0.411(36) 


-0.175(36)(17) 


d 


-0.203(35) 


0.205(35)(0) 



at next-to-leading order in the MS(MS) scheme. 

An overview of lattice QCD calculations of nucleon 
spin observables, in particular of moments of generalized 
parton distributions that give access to the total quark an- 
gular momentum Jq, can be found in |34| . Here we focus 
on the latest published results from the LHP collabora- 
tion |35) . They were obtained in the framework of a mixed 
action approach with np = 2-1-1 dynamical fermions, with 
lattice pion masses as low as « 300 MeV. The compu- 
tationally demanding quark line disconnected diagrams, 
which contribute in the singlet sector, were not included 
in this study. The final values for ASq, Lq and Jq at the 
physical pion mass were obtained from extrapolations em- 
ploying the covariant baryon chiral perturbation theory 
results of |36j . We refer to the original publication 35 for 
the details of the lattice simulation, the numerical analy- 
sis, and a discussion of the statistical and systematic un- 
certainties. A summary of the lattice results, for the MS 
scheme at a scale of 4 GeV , is given in Table 2 The errors 
given in this table do not include systematic uncertainties 
from chiral extrapolations and disconnected diagrams. 

For our QCD evolution of lattice results, we assume a 
vanishing contribution from strange quarks (see also |37| ). 
The total gluon angular momentum is given by Jg = 
\ — Jri+d- As starting values at /i^ = 4GeV^, we use the 
numbers given in Table [2] For the running coupling, we 
set a^^^^^''^='^{Mz) = 0_A184 |38j, and employ the fia- 

NNLO,nF=3 _ 



vor matching conditions 



0.338 GeV, A 



NLO,nF=3 



39 



40 



^QCD 

GeV. With this input 



to obtain ^qcq 
(1388GeV, and A^^^''=^ 



0.148 



we have solved the LO, NLO and 
NNLO coupled evolution equations and found the scale 
dependence plotted in Figures [T] and [2] 

The results at LO, NLO and NNLO, employing the 
standard analytical expressions for the perturbative strong 
coupling constant (corresponding to an expansion in 
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Figure 1. Scale dependence of AS/2, Lq, Jg, and Lu-d shown together with ^qa = ^AEu-d, starting from the lattice QCD 
results at /i^ = 4 GeV^ given in Table [2] In all diagrams the solid, dashed and short-dashed black lines are solutions of the QCD 
evolution equations at NNLO, NLO, and LO, respectively. The colored bands are obtained by imposing upper bounds for Q3 
(see text). 
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Figure 2. Evolution of L„ and Ld in NNLO (solid lines), NLO 
(dashed lines) and LO (short dashed lines). The shaded areas 
are assigned as in Fig. [l] 



ll\n{y?/A^) beyond LO, see, e.g., jsSj) in the MS scheme 
at the appropriate order, are given by the short-dashed, 
dashed, and solid black curves, respectively. We note that 
the deviation of the approximate analytical expressions 
for as from the exact (numerical) solutions of the evolu- 
tion equations increases as one approaches lower scales. 

For rii. = 3 and ^qcd°'"''"^ = 0.338 GeV, the formally 
exact solution for the running coupling at NNLO would 
already diverge around /i^ '-^ 0.27 GeV^. The curves in 
Fig. [1] obtained for as in the 1/ ln(/i^/yl^)-approximation 
are therefore only indicative for a strong coupling constant 
that grows indefinitely as /i^ — > 0. 

A comparison with the model results, e.g. as proposed 
by Myhrer and Thomas [7 - 10 , requires an evolution down 
to scales /i^ ^ 0.1 — 0.3 GeV , far away from the pertur- 
bative QCD regime. From the results in Figs. [T| and [2j we 
find that the evolution curves at NLO and NNLO begin 
to show a very strong curvature exactly in this region. 
Clearly, at such low scales quantitative statements based 
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b.) 





d.) 




Figure 3. Qs-dependence of the momentum and angular momentum contributions {x)u-d and Lu-d with starting values fixed 
at Usfi = 0.28. Diagrams c and d show the corresponding renormalization scale dependence. In all the diagrams solid, dashed 
and short-dashed red lines correspond to NNLO, NLO, and LO, respectively. 



on a perturbative QCD analysis (including the running of 
as) are no longer reliable. 

With respect to (the non-perturbative) ag, one would 
expect in any case that it saturates at low scales, as sug- 
gested by non-perturbative resummation in the infrared 
region |41j-|43j. A further rough impression about the un- 
certainties in the evolution may therefore be obtained as 
follows. As an alternative to the infrared divergent, per- 
turbative coupling UsifM^) in the evolution equations, we 
consider an effective a®*^(/x) that approaches a fixed value 
^s'max small /i^. For the corresponding numerical cal- 
culation we use af^{^^) — as{is^) of appropriate order 
in the MS scheme for all /i for which as(/i^) ^ Oiflnax- 
Below the scale fig at which Q!s(/io) — as^iax; ^^^^ 
af^ = Qfg'^ax. For illustration, we chose two different val- 
ues, a^^ax — 1-5 8i,nd a^^ax — 3. Performing the evolu- 
tion with these restricted couplings spans the shaded col- 
ored bands in Figures [T] and [2] The boundary with flatter 
evolution always corresponds to a^^jax = 1-5, the steeper 
one to a^^jax = 3- The lighter colored bands correspond 
to NLO, the darker colored bands to NNLO. The rela- 
tively small j^qq'j^^''"'^ = 0.148 GeV, obtained from the 
flavor matching procedure, implies that the correspond- 



ing Q!s(/i^) at LO stays below Q;°*[„ax ™ the considered 
region of /i^. Hence no bands are shown in this case, and 
one flnds the remarkably stable (but unrealistic) LO evo- 
lution shown in the flgures. 

NLO and NNLO evolution results in the dashed and 
solid lines, and the lighter and darker shaded colored bands, 
respectively. As mentioned before, in contrast to the LO 
evolution, strong evolution effects can already be seen at 
scales fj? < 0.3 GeV^. Before discussing the results in more 
detail, we note as a general feature that the bands ob- 
tained for our choices of a^^ax start to broaden quickly 
below ~ 0.3 GeV^, indicating potentially large uncer- 
tainties in the evolution as one enters the non-perturbative 
regime. With the exception of AS, we observe a broad 
overlap of the NLO and NNLO bands for each of the dif- 
ferent observables. 

At these orders we find interesting qualitative and quan- 
titative changes of the proton spin decomposition under 
evolution. The singlet quark spin contribution AS be- 
comes scale dependent at NLO and increases with lower 
scales, a behavior that is even more strongly pronounced 
in NNLO. Similarly, the contribution from Lg stays posi- 
tive and grows quickly at low scales, while Jg crosses zero 
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in the region of ^2 ^ 0.21 to 0.29 GeV^ and then becomes 
large and negative. At the same time, shows a 

strong upwards bending and moves from its negative start- 
ing value at higher scales towards large positive values at 
low scales, crossing zero around /i^ ~ 0.28 to 0.32 GeV^. 
From the evolution of Lg and L^-d one can deduce the 
separate /i-dependences of L„ and L^, as displayed in Fig- 
ure [2j Both Lu and Ld change sign under evolution. As 
one moves in the direction of lower scales, the contribution 
from the up-quarks, L„, changes from negative to increas- 
ingly large positive values at about /i^ ~ 0.32 — 0.35 GeV^, 
while the zero crossing of from positive to increasingly 
large negative values takes places at slightly lower scales 
of fj,^ ^ 0.2 to 0.3 GeV^. Similar trends have already been 
observed in the LO-study of Ref. [s] , based however on a 
substantially larger ^qcd^~'^ = 0.24 GeV. It is interesting 
to observe that the crossing points of L„ and Ld roughly 
match as we proceed from NLO to NNLO. 

For a better understanding of the scale dependence at 
very low scales, i.e. for large values of the strong coupling 
constant, it is instructive to present the evolution in terms 
of as instead of ^j^This is illustrated in Figure [3] for the 
nucleon isovector momentum fraction {x)u-d on the left 
and for L^^d on the right. Note that the evolution equa- 
tions for both quantities are based on the same anomalous 
dimension. The results were obtained by rewriting the evo- 
lution equations, i.e. translating derivatives with respect 
to fj? into derivatives with respect to as, and then re- 
placing ^^°^2 by the QCD /3-function. As starting values 
we have used {x)u-d — 0.155 and L^-d = —0.38 [35^ at 
asfi = 0.28 (corresponding to a scale of /i^ ^ 4GeV^). 
Both quantities are remarkably stable under evolution in 
as, even at large as - For example, at as = 2.0 the momen- 
tum fraction {x)u-d at NLO is about ~ 14% larger than 
the LO result, and NNLO and NLO results differ by only 
about ~ 8%. Apart from the difference in the starting val- 
ues, the form of the as-dependence is identical for L^-d 
and {x)u-d- The zero crossing of Ly^^d quickly moves to 
smaller as as the order of the perturbative evolution is 
increased. Figures [Sj; and d show how the as dependences 
translate to a renormalization scale dependence. The log- 
arithmic dependence of as{n) that shows up strongly as 
fi approaches /Iqcd , produces the strong curvature in the 
evolution of the observables as functions of the scale fi. 

Small but unimportant differences between Fig. 3d and 
Fig. Id result from the slightly different procedures in- 
volved, as described. 

Further remarks on the evolution and a comparison 
of evolved lattice results with calculations performed in a 
chiral quark model will be presented below in Section |4] 



^ showing less systematic uncertainties in the lattice compu- 
tation as contributions from disconnected diagrams cancel out 
for isovector quantities 

^ Ph.H. would like to thank M. Diehl for helpful discussions 
on this point. 



3 Contributions to the nucleon spin in a 
chiral quark model 

3.1 Pion cloud contributions, revisited 

It is a well established fact that the nucleon is a com- 
plex many-body system, with the three valence quarks 
and multiple quark-antiquark pairs embedded in a strong, 
non-perturbative gluonic field configuration. Chiral quark 
models draw a simplified picture of this complexity in 
terms of valence quarks in a confining bag coupled to the 
pion cloud, based on spontaneously broken chiral symme- 
try in low-energy QCD. A frequently used representative 



of such chiral models is the cloudy bag (CBM) [44| 
that couples the pion cloud to quarks in the MIT bag 
such that chiral invariance is realized in the limit of mass- 
less quarks. This section summarizes the present status 
concerning nucleon spin structure from this model point 
of view. 

The relativistic treatment of quarks itself yields al- 
ready results that differ significantly from the ordinary 
SU(6) quark model predictions. The AS = 1 of the non- 
relativistic quark model is reduced to about AE^^'^ = 
0.65 in the MIT bag model. The "missing spin" is in- 
terpreted as orbital angular momentum of the valence 
quarks, 2L^J^ — 0.35, associated with the lower compo- 
nents of the Dirac quark wave functions. 

The correction factors for the pion cloud in the CBM 
were already derived by Myhrer and Thomas in j7 - 10 
For the singlet expectation values. 



AS 



CBM 



0.65 • 7Ts(i?), 2L^™ = 0.35 • 77s (i?), 



27. 



CBM 



= l-7Ts(7?), 



(12) 



and for the non-singlet expectation values 37 



zii:„c_T=g(f)^«^ = |.0.65 . 77^s(i?), 



27.: 



CBM 

i—d 



0.35-7Tns(7?). 



(13) 



We have denoted the pion cloud correction factors by 
77s (7?) and 7Tns(7?), each for a given bag radius 7?. For 
their explicit form we refer to |7 ■ 10 37 and references 
therein. We have reproduced these factors using the for- 
malism described in f46' . Their radius dependence is plot- 
ted in Figure [4] The singlet correction factor, 77s, is smaller 
than unity and hence leads to the expected reduction of 
the quark spin contribution. At the same time, 7Tns < 1 
leads to a slightly less favourable comparison of g^-* 
with the experimental value of qa- This mismatch is a 
feature that depends on the chiral representation cho- 



sen for the model. In particular, as pointed out in 37 



choosing a volume coupling version instead of the stan- 
dard surface coupling reproduces the experimental value 
of with very good accuracy. A common way of treat- 
ing the discrepancy between g^^'^^^^ and the empirical 



9A = 9 a 



(3) 



1.27 is by inclusion of a phenomenological 
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Figure 4. Radius dependence of the singlet and non-singlet 
corrections 77s and 77ns associated with the pion cloud of the 
nucleon. 



center-of-mass correction. This correction is just a mul- 
tiplicative factor for AS and Z\Z't(_(;. For consistency, 
however, one should then also rescale Lq and Lu_ii ac- 
cordingly to keep the spin sum rule conserved. The cor- 
responding results are shown in Table [3] in addition to 
results obtained without the center-of-mass corrections. If 
explicit gluon ope rators are taken into account (see Table 
Inland section 3.2 below), these corrections cannot be ap- 
plied uniquely. In that case we give only the model results 
without rescaling by cm. corrections. 

A detailed analysis of AS also requires, in principle, 
a discussion of g\ since the flavour singlet axial vector 
coupling constant extracted from polarised deep inelas- 
tic scattering is sensitive to that value. In the present 
case we restrict ourselves to a flavor-SU(2) cloudy bag 

model which implies AS — g\ . This is quite compatible 
with the small strange quark contribution As ^ —0.01 
discussed in Ref. iSTi. 



3.2 Corrections from one-gluon exchange processes 

The MIT bag model produces degenerate masses of the nu- 
cleon and the delta, whereas the empirical mass splitting 
is about 300 MeV. In order to account for this mass differ- 
ence, an additional spin-spin interaction between quarks 
is introduced. The traditional ad-hoc way of doing this 



Notice that the value = 1.27, given in [7 



10 



is ob- 
tained by adjusting the phenomenological center-of-mass cor- 
rection, while this correction has not been included for any of 
the other spin observables listed in these references. 
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Figure 5. One-gluon exchange (OGE) corrections for AS, 
Ls,NS and g^^^ ■ Diagram (a) and (b) are contribution from in- 
termediate quark states and (c) and (d) from quark-antiquark 
pairs. 



is by allowing one-gluon exchanges between quarks in the 
interior of the bag, with an effective coupling fig. This cou- 
pling should not be confused with the strong coupling of 
QCD. It represents a free parameter chosen in such a way 



that the model reproduces the light hadron spectrum 48 
The values used in the literature vary between as = 1^491 
and as — 2.2 \50 . We use these two values as options in 
the results shown later in Table |3] and |4l 

We have performed calculations in analogy to Ref. [50] , 
where the color magnetic corrections to baryon magnetic 
moments and to semileptonic decays, i.e. the axial cou- 
pling constant, were derived at order cts- Following the 
arguments given there, we neglect the color electric cor- 
rections and drop loop diagrams. That means, we consider 
diagrams [5f^a)-[5]^d) , in which only color magnetic gluons 
are exchanged. 

For the singlet expectation values, AS and Lq, we find 
the (additive) OGE corrections 



SAS = -2dg ■ as 



SLq = Sg ■ as 



(14) 



with (5g ~ 2.5 • 10 ^, where Lq is used in its non-gauge- 



invariant formulation (|3|, and intermediate (anti-)quarks 
in the orbitals P1/2 , F3/2 , £'3/2 , > ^'i'/2 > ^3/2 ' -^3/2 ' ^2 
are taken into account (conventions are chosen as in [50]). 
As already pointed out in [7}|l0] the corrections are mainly 
due to antiquarks propagating in the /'i/2,^3/2 orbitals. 
Compared to SAS_ 0.15 and SLg - 0.08 for 5^ = 2.2 



as presented in |8j[9], our corrections are slightly smaller. 
For the non-singlet operators we find: 



6ASu. 



X (3) 
^9 A 



as. 



(15) 
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At this point our model agrees with that of Refs 
For comparison, explicit numbers for the singlet am 
singlet contributions to the nucleon spin in the MIT bag 
model, as well as the OGE-improved MIT bag and cloudy 
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Table 3. Spin structure of the nucleon in the MIT bag model, with corrections from one-gluon exchanges (OGE), from the pion 
cloud and from center-of-mass rescaling. The non-gauge- invariant decomposition of the nucleon spin, Eqs. ([3|, is used here. 





AE/2 


L, 


AE^^a/2 




relativistic (MIT bag model) 


0.33 


0.17 


0.54 


0.29 


+OGE (qs = 1.0): 


0.30 


0.20 


0.55 


0.28 


(5s=2.2): 


0.27 


0.23 


0.56 


0.27 


+pion cloud (R — Ifm, as = 1.0): 


0.24 


0.26 


0.51 


0.26 


(R = Ifm, as=2.2): 


0.22 


0.28 


0.52 


0.25 


+center of mass {R = Ifm, as=1.0): 
{R = Ifm, &s=2.2): 


0.30 
0.27 


0.20 
0.23 


0.64 
0.64 


0.20 
0.21 



Table 4. Spin structure of the nucleon, based on manifestly gauge invariant operators, in the MIT bag model, together with 
corrections from one gluon exchanges (OGE) and from the pion cloud. 





AE/2 


rGI 


■^s 


AEu-d/2 


rGI 

^u~d 


relativistic (MIT bag model) 


0.33 


0.17 





0.54 


0.29 


+OGE (5, = 1.0): 


0.30 


0.40 


-0.20 


0.55 


0.21 


(a,=2.2): 


0.27 


0.68 


-0.45 


0.56 


0.12 


+pion cloud {R = Ifm, as = 1.0): 


0.24 


0.42 


-0.16 


0.51 


0.19 


(i?= Ifm, &s=2.2): 


0.22 


0.64 


-0.36 


0.52 


0.10 



Figure 6. One-gluon exchange diagram attached to a quark- 
quark-gluon interaction vertex. 



Furthermore, from the gauge- invariant spin sum rule Eq. (|4| 
we conclude that the contribution from the total gluon an- 
gular momentum equals 



J?' = -SL, 



q,A- 



(18) 



bag model, for two different values of as, are displayed in 
Table |3l0 



Once explicit gluon degrees of freedom, e.g. in form 
of one gluon exchange processes, are taken into account, 
questions of gauge invariance of the calculation must be 
carefully addressed. For a consistent calculation that in- 
cludes gluon exchange contributions at order dig, and for a 
meaningful comparison with results from lattice QCD, we 
have to employ the gauge-invariant orbital angular mo- 
mentum operator L^^ in Eq. ^ instead of Lq defined 
in (jsj). The covariant derivative in L^^ produces an addi- 
tional quark-gluon interaction so that one must take into 
account the diagram in Fig. [6] for the corrections at or- 
der cts to the MIT bag expectation values. This diagram 
yields the large contribution 



q,A 



0.203 5^ 



(16) 



where the subscript A stands for the gauge field interac- 
tion term. The total correction to the quark orbital mo- 
mentum is then given by = SLq+SLq A- The diagram 



in Fig. 

shifts i 



also contributes to the gauge invariant L^^^ and 



5Lu^ 



d,A 



= -T^SLq^A- 



(17) 



* We tabulate these results here only for historical reasons 
as they are based on the non-gauge-invariant decomposition, 



Eq. ^. 



We notice that the corrections (16l-(18l are much larger 
than the known one-gluon exchange contributions from 
the diagrams of Fig. [S] given in Eqs. ( [T4| , (15 1. In par- 
ticular, SLq A dominates L^'^"^ for the chosen parameters. 
The MIT bag model results for the gauge invariant decom- 
position of the nucleon spin are summarized in Table [4] 
together with the combined results, including relativistic 
effects plus one-gluon exchange corrections plus correc- 
tions from the pion cloud. 



To conclude this section, we emphasize that a direct 
calculation of J^^ = {P + \ J d^r[r x (E x B)]^\P+) (i.e. 
not invoking the spin sum rule) in the framework of the 
model requires a careful treatment of the boundary condi- 
tions for the color electric fields. The boundary conditions 
f • E\r=R cannot be fulfilled for the color electric fields, in 
the way described in 51 . This leads to a non-vanishing 
surface term in the calculation of J^^ and, therefore, po- 
tentially to inconsistencies with respect to the spin sum 
rule. A calculation with color electric fields that obey the 
boundary conditions, as given in |52| , turns out to be sig- 
nificantly more involved and will not be described in this 
work. Note, however, that SLq^A is not affected by such 
complications since the corresponding operator does not 
involve color electric fields. It is therefore legitimate to 
extract the corresponding gluon angular momentum from 



When J^^ is calculated directly, using the "wrong" 
color electric fields, it spoils the spin sum rule. Actually 
the direct evaluation of J^^ can be used to check our result 



M. Altenbuchinger et al.: Spin structure of the nucleon: QCD evolution, lattice results and models 



9 



for SLq A- Consider the decomposition 

J (fxx X {E X B) = j d^xEx A + j d^xE'{x x \/)A' 

-J Sxgip\x X A)^l}- J d^x\7'[E^{x x A}]. (19) 

The left hand side corresponds to J^^, the right hand 
side to AG + Lg — 5Lq A supplemented by a surface term, 
— / d?xV^[E^{x X A)]. This surface terms vanishes in the 
free field theory but in our model calculation this is not the 
case. Therefore, the total gluon angular momentum calcu- 
lated through the spin sum rule equals J d^xx X [E X B) + 
J d'^x\/^[E^ {x X A)], which indeed can be confirmed by 
a direct calculation. In fact, the OGE corrections to AG 
and Lg cancel each other. 

4 Discussion and summary 

The present study has been motivated by the observation 
of an apparent contradiction between quark orbital an- 
gular momentum contributions, Lg and Ly^_ii, calculated 
in models and derived from lattice QCD computations. 
At the same time, model and lattice QCD results for the 
quark spin contributions AS and = ASy^^^ are rea- 
sonably consistent once pion cloud and gluon exchange 
effects are incorporated in the model |7|. When compar- 
ing the two approaches, it is essentiaT to note that all 
spin observables, except gA, are scale (and scheme) de- 
pendent quantities. Since the model scales are typically 
low, ~ 0.1 — 0.3 GeV^, a careful study of the scale evo- 
lution is necessary. In contrast to previous studies ^ , our 
investigations are based on a "downwards" evolution of the 
lattice results, starting at higher scales where a perturba- 
tive treatment appears safe, to low (model) scales where 
higher-order effects must be taken into accout. By consid- 
ering the evolution at LO, NLO and NNLO, together with 
the possibility that the (non-perturbative) strong coupling 
saturates at very low scales, this approach allows us to 
perform a meaningful comparison of the lattice results 
and the spin contributions obtained in different model ap- 
proaches. 

With the comparatively low value of A^^''='^ = 0.148 
GeV obtained from the flavor matching procedure, the LO 
evolution is rather flat for all observables down to scales 
of /i-^ = 0.1 GeV^. At this stage it is not possible to resolve 
the aforementioned contradiction between phenomenolog- 
ical or lattice results and the model results from Table [H 
This is in contrast to the observation in Ref. |8j. 

Extending the evolution equations to NLO and NNLO, 
one flnds a strong /^^-dependence at typical model scales. 
Even if the matching scale /^^odci small to draw 

quantitative conclusions, the trends are indicative. Em- 
ploying different values for ^qcd'"''"'^ ^^'^ ^qcd'^'"''"'^ 
as obtained from flavor matching, the NLO and NNLO 
results for Jg, as well as for L^-d, show very good over- 
lap even at the lowest scales. For AS and Lq, the results 
at NLO and NNLO are still quantitatively comparable in 



magnitude down to scales of about /i^ ^ 0.3 GeV , but 
show larger deviations as /i — O.lGeV^ In the region of 
such low scales where the evolution effects become strong, 
there is indeed an overlap with typical model results for 
the individual observables, as given in Table |3] and |4] 

Incidentally, the NLO evolved lattice results for AS, 
Lq and Lu-d turn out to be reasonably close to the orig- 
inal MIT bag values at a scale of ^i^adc\ ^ 0.28 GeV^ 
where the total contribution from the gluons vanishes, 
JgiA,odc\) = 1/2 - Lq - AS/2 ~ 0. As the model calcu- 
lations are improved, however, this apparent consistency 
deteriorates. Inclusion of the pion cloud in the CBM low- 
ers AS and increases Lq signiflcantly. The inclusion of 
the phenomenologically center of mass correction, which 
is not based on solid theory, would in part restore the 
agreement. On the other hand, inclusion of further gluon 
exchange corrections as in Table [4] would make the match- 
ing with the evolved lattice data at a common low scale 
progressively more difficult. For example, while model im- 
provements lead to a significantly larger L^^ , thereby im- 
plying a lower matching scale, the corresponding results 
for the isovector L'^^^ become successively smaller, which 
would in turn require a matching at an increasingly higher 
scale. 

It might seem that the scale dependence of AS im- 
plies a large gluon spin fraction AG at the scales of po- 
larized DIS, larger in magnitude than the constraints pro- 
vided by the HERMES, COMPASS and RHIC measure- 
ments [5 20 - 22 . As a test we performed the NLO evolu- 
tion downward starting from AG — Q &i — 4GeV^. 
One then finds AG ^ —0.4 around the scale where Jg 
vanishes. 

Concerning the lattice calculations, one source of sys- 
tematic uncertainty can be eliminated by studying isovec- 
tor quantities such as Ly_d for which disconnected dia- 



grams, not taken into account in Ref. (35 
From the model investigations (see Tables 



pects in the range of 0.1 ^ 

lattice results start negative at fi^ ^ 4 GeV^ . We find that 
the downward QCD evolution does indeed predict the ap- 
propriate change of sign (see Fig. [iji) at NLO and NNLO. 
Other systematic uncertainties on the lattice side, for ex- 
ample those related to lattice operator renormalization is- 
sues, would affect the normalization of L^-d but would not 
change this picture significantly. As already noted above, 
and in contrast to the singlet Lq, inclusion of explicit gluon 
degrees of freedom in the properly gauge invariant treat- 
ment of the quark orbital momentum operator leads to a 
reduction of L„_rf at model scales (see Tables [3] |4]) and 
moves this quantity closer to the extrapolated lattice QCD 
results. The sign change of Lu~d can be traced in detail by 
examining the crossing of L„ and Ld as shown in Fig. [2] 

In summary, our analysis underlines the difficulty of a 
simultaneous, quantitative understanding of model calcu- 
lations and lattice QCD results for the decomposition of 



0.3 



a 



cancel out. 
3 1 |4l one ex- 

the 



In contrast 

2 



^ Note, however, that a recent calculation using a chiral 
quark soliton model gives a negative L^-d even at low scales 



53 . 
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the nucleoli's spin into the angular momenta of the con- 
stituents. While the perturbative corrections from NLO 
to NNLO for the evolved lattice results are at a tolerable 
level, the broad bands we obtain for a non-perturbative, 
saturated as indicate potentially large systematic uncer- 
tainties in the evolution at very low scales. On the side of 
the model calculation, we find that the effects from dif- 
ferent types of improvements (related to the pion cloud 
effects, manifestly gauge invariant 0AM operators, one- 
gluon exchange corrections) tend to make it increasingly 
difficult to find a common low matching scale where at 
least a semi-quantitative agreement with the evolved lat- 
tice results can be achieved for all of the different spin ob- 
servables. Conversely, this indicates that it will be difficult 
to arrive at quantitatively reliable predictions from model 
calculations starting at scales smaller than ^ 0.3 GeV^ 
and evolving upward to scales accessible in experiments 
and related QCD phenomenology. We stress that this ob- 
servation is not in contradiction with the phenomenologi- 
cally very successful approach of "dynamically generated" 
(unpolarized) parton distributions (PDF) |54-56 , where 



specific ansatze for the x-dependent PDFs are evolved 
from a very low initial scale, /j,o < 1 GeV, to higher scales. 
The latter approach necessarily involves adjustable pa- 
rameters (for each type of PDF) to achieve a fully quanti- 
tative describtion of the experimental DIS data at scales 

= /x^ > 1 GeV^ and over a wide range of the momen- 
tum fractions x. 

A possible exception concerning the previous critical 
assessment is the isovector orbital angular momentum com- 
bination Lu-d for which systematic lattice errors are min- 
imal. This quantity displays generic behavior with a sign 
change as it evolves from lattice QCD to low scales, in 
accordance with the model expectations listed in Tables 
[3] and |4j The stability of this evolution, as one proceeds 
from NLO to NNLO, becomes particularly apparent when 
plotted as a function of and compared with the cor- 
responding evolution of the average momentum fraction. 
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